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Nuclear spins in the solid state environment of diamond are highly coherent, but difficult to 
rapidly control due to the small nuclear gyromagnetic ratio. Here we demonstrate a more than 
50-fold enhancement of the effective nuclear gyromagnetic ratio by coupling the nuclear spin to an 
electronic spin of a nitrogen-vacancy (NV) center in diamond. The enhancement allows for faster 
nuclear spin rotations and is in good agreement with second-order perturbation theory. The method 
may be applied to other systems with similar electron-nuclear spin interactions, such as phosphorous 
donors in silicon, opening up the possibility of fast and direct nuclear spin control in coupled spin 
systems. 

PACS numbers: 03.67.Lx, 31.30.Gs, 76.30.Mi, 76.60.-k 


Advancements in quantum computation have been 
accelerated by the development of qubits with long 
quantum coherence times and high fidelity quantum 
control^. The nuclear spin degree of freedom has long 
been considered to be a good candidate for a qubifcA— . 
Several approaches towards solid state quantum compu¬ 
tation rely on using the nuclear spin as a memory and 
impressive coherence times have been achieved^—. How¬ 
ever, many of these approaches rely on using an auxiliary 
qubit for manipulation of the quantum state, as direct 
control of the nuclear spin is typically limited to 10-50 
H s Rabi periods due to the small nuclear gyromagnetic 
ratio 11 ! 12 . In comparison, electron spin Rabi periods as 
short as 4 ns have been achieved in a variety of condensed 
matter systems^ - — . 

Recent experiments on nitrogen-vacancy (NV) centers 
have shown that nuclear spins coupled to a NV center 
can exhibit Rabi oscillations that are significantly faster 
than expected for a bare nuclear spin that is coupled 
to an ac magnetic fiel d 17 ! 18 . The faster rotation rate 
can be interpreted as an enhancement of the effective 
nuclear gyromagnetic ratio that results from hyperfine 
interactions with the NV electronic spin. 

In this Rapid Communication, we show that the en¬ 
hancement of the nuclear spin gyromagnetic ratio is tun¬ 
able by varying the electronic spin transition frequency. 
We demonstrate over a factor of 50 improvement in the 
nuclear spin rotation rate compared to what is expected 
from the bare nuclear gyromagnetic ratio^. The en¬ 
hancement is observed by directly measuring the nu¬ 
clear spin Rabi frequency and is in good agreement with 
second-order perturbation theory^. The enhanced nu¬ 
clear spin Rabi frequency may allow nuclear spins to be 
utilized in quantum control gate sequences, in addition 
to being a resource for a long-lived quantum memory. 

The quantum system consists of a NV electron spin 
and its intrinsic 14 N 1 = 1 nuclear spin. A naturally 
occurring NV center is selected from a liigh-purity type- 
Ila diamond sample (Element Six). The experimental 
setup is schematically illustrated in Fig. Oja). A suit¬ 
able NV center is located using a confocal microscope 
with a 532 nm excitation laser. Photoluminescence (PL) 



FIG. 1: Experimental setup: (a) A confocal microscope is 
used to access a single naturally occurring NV center in high- 
purity diamond. The NV electronic spin is optically initialized 
and read out with a 532 nm laser, while the intrinsic 14 N of 
the NV is used as the nuclear spin qubit. Quantum control is 
achieved by applying MW and RF pulses through the short- 
terminated stripline fabricated on the surface of the diamond, 
(b) Circuit diagram: Two signal generators are used to gen¬ 
erate MW and RF pulses. The signals are combined with an 
adder before being delivered to the stripline. 

from the NV center is collected using a high numerical 
aperture (NA = 0.95) objective and directed towards sin¬ 
gle photon detectors using a combination of fiber and 
free space optics^. An ac magnetic field with amplitude 
B ac is generated using a stripline that is fabricated on 
the surface of the diamond. Microwave (MW) and ra¬ 
dio frequency (RF) pulses are applied to the stripline to 
drive electronic and nuclear spin rotations. The circuit 
diagram is shown in Fig- HKb) . A SRS SG394 (Aglient 
N5181A) signal generator is used to generate MW (RF) 
pulses. The MW signal is amplified with a broadband 
amplifier (Triad RF TA1003) to allow for fast manipu¬ 
lation of the electronic spin. RF and MW signals are 
combined with a resistive splitter-combiner before they 
are delivered to the stripline. By driving electron Rabi 
oscillations and measuring the resulting Rabi frequency, 
we estimate I?r,f ~ 2.9 G at the sample with 1 mW of 
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FIG. 2: (a) NV center energy level diagram. A selective MW 
pulse is tuned to resonance with the | ms, mi) = |0, +1) to | — 
1,+1) transition for initialization and readout of the nuclear 
spin, (b) PL as a function of <5mw = i'MW — A, showing 
nearly full nuclear spin polarization in the mj = +1 level at 
the ESLAC (lower trace, B z = 514 G). In comparison, the 
polarization is significantly reduced away from the ESLAC 
(upper trace, B z = 847 G). The detuning <5mw is defined 
relative to the center frequencies 1429.7 MHz (496.4 MHz) 
for B z = 514 G ( B z = 847 G). Data sets are vertically offset 
for clarity, (c) Pulse sequence used for nuclear spin Rabi 
oscillations (see text for details), (d) PL as a function of RF 
pulse duration trf showing nuclear spin Rabi oscillations. 


power applied to the stripline. We verify the accuracy of 
this conversion across the range of frequencies and pow¬ 
ers used in this experiment by performing an independent 
multi-frequency excitation^ on the electronic spin (data 
not shown). By simultaneously driving the NV with RF 
and MW pulses, the slowly varying RF field periodically 
shifts the NV electronic levels, resulting in a broadened 
MW transition. From the width of the MW transition, 
we can extract the value of £?rf using the well-established 
electron g-factor— g e = 2. 

The NV center Hamiltonian H = Hq + V results in the 
level diagram shown in Fig. 2(a) and consists of secular 
terms Hq, and non-secular terms V: 


H 0 /h = DS Z + g e ^BB z S z - gNl^wB ■ I 

— S z A Z jIj + QI z , ( 1 ) 

j=x,y,z 

V/h = i (g e fi B (B_S + + B+S-) 

- ^2 {S+A-jlj + S-A+jlj) 

j=x,y,z 

( 2 ) 

Here g e ( gw ) is the electronic (nuclear) g-factor, y B (hn) 
is the Bohr (nuclear) magneton, B is the external mag¬ 
netic field, and S (I) are the electron (nuclear) spin oper¬ 


ators. Hyperfine couplings are given by the matrix com¬ 
ponents Aij , where i,j = {x,y : z}. Raising and lowering 
operators are defined as B± = B x ± iB y , S± = S x ± iS y , 
A±j = A x j ± i A y j . We neglect interactions with nearby 
13 C impurities as their interaction strengths are much 
smaller than the coupling to the intrinsic 14 N and no 
such coupling is visible in the spectroscopy data [Fig. 
[21(b)] . Electronic spin states ms = 0 and ms = —1 are 
separated by A — D — g e fi B B z , where D = 2.87 GHz 
is the ground state zero-field splitting and g e fi B B z is the 
Zeeman shift from the external magnetic field B z applied 
along the NV symmetry axis (g e yB = 2.802 MHz/G). 
The ms = +1 spin state is far detuned and therefore not 
shown in Fig. [21(a). For each electronic spin state, the 
energy levels are further split into three sublevels by the 
14 N quadrupole coupling Q = —4.96 MHz and nuclear 
Zeeman shift (gNBN = 0.308 kHz/G). Electron-nuclear 
axial hyperfine coupling A zz = —2.16 MHz further sep¬ 
arates the mi = —1,+1 states in the ms = —1 sub¬ 
space. We note that while the perpendicular hyperfine 
terms A± = A xx = A yy = —2.70 MHz are non-zero, 
the electron-nuclear flip-flop process is suppressed due to 
large quadrupole coupling Q. The off-diagonal terms of 
the matrix A are zero due to the symmetry of the NV 
cente r 9 ' 17 ! 21 . 

We select the two well-isolated sublevels | ms,mi) = 

| — 1,+1), | — 1,0) for the nuclear spin qubit to al¬ 
low for selective excitation using the RF field (frequency 
i/ RF ~ 3 MHz). This choice of states for our qubit 
also provides a convenient means of nuclear spin read¬ 
out. We can directly map the nuclear spin state to the 
electronic spin state by applying a selective microwave 
7r-pulse (with frequency ^mw tuned to resonance with 
To/ = +1 transition) at the end of the nuclear control 
sequence. The electronic spin state is then read out by 
optical excitation^ill. 

For the experiment, B z is aligned to within ±1 degree 
of the NV-axis22. We choose to work close to the excited 
state level anti-crossing (ESLAC), B z ~ 500 G, to allow 
for efficient polarization of the nuclear spin to the mi = 
+1 state by optical pumping with a 532 nm lase r 23 ! 24 . To 
probe the electronic spin transition frequencies, we ap¬ 
ply pulsed microwave excitation with varying frequency 
l'mw ■ When z/mw is on resonance with an electronic tran¬ 
sition we observe a dip in the PL intensity. Figure [21(b) 
shows the PL as a function of z^mW; showing the tran¬ 
sition frequencies associated with different nuclear spin 
projections^. Two data sets are shown for comparison. 
At the ESLAC ( B z = 514 G), the nuclear spin is fully po¬ 
larized to ms = +1 (polarization P > 95%). Away from 
the ESLAC the polarization is less efficient, as shown in 
the data with B z = 847 G (P ~ 72%), where transitions 
corresponding to mi = —1,0 can also be observed. While 
imperfect polarization reduces the contrast of the nuclear 
spin Rabi oscillations, it does not affect the measurement 
of the Rabi frequency and the results that follow. 

The full experimental pulse sequence is illustrated in 
Fig- Oc). We start by optically pumping the electron 
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FIG. 3: Nuclear spin Rabi frequency fljv plotted as a function 
of S RF for several different values of B z . The slope of each 
line corresponds to the effective nuclear gyromagnetic ratio 
7jv,eff- For the largest value of B z = 873 G, 7iv,eff/7iv = 57. 


by Eq. |T] and Eq. [2J respectivel y 19 ! 26 . Treating V as a 
perturbation, and calculating to second order in pertur¬ 
bation theory, we have: 



(3m 2 s -2)D-m s geHBB z 

2(D*-g*> l %B*) [M+ " h 


( 3 ) 


where 


M = -g e Li B (B-A + ■ 1 + B+A _ • /), (4) 


and 


N = g 2 e n 2 B B 2 ± + 2A + A_+(A+ x !_)•/. (5) 


Here Aj = (Aj X , Aj y , Aj Z ) and A± = A x ± iA y . We can 
neglect all the terms in N as they commute with I x and 
I y and therefore do not contribute to spin flips. The only 
remaining term, given in Eq. [G] contributes directly to 
the nuclear spin Rabi frequency and the effective nuclear 
gyromagnetic ratio. For ms = — 1, we have, 


and nuclear spins with a 4 -/is long 532 nm laser pulse. 
After the system is polarized into the |0, +1} state, a 
selective MW 7r-pulse is applied to transfer the popula¬ 
tion to the | — 1,+1) state, completing the initialization 
process. We then drive nuclear spin Rabi oscillations by 
applying a RF pulse with varying duration trf resonant 
with the | — 1,+1) to | — 1,0) transition. Finally, optical 
readout is performed by applying another selective MW 
7r-pulse that converts the population from | — 1,+1) to 
the bright state |0, +1). This yields a PL signal that is 
proportional to the | — 1,+1) population at the end of 
the pulse sequence. 

Nuclear Rabi oscillations are readily observed by im¬ 
plementing this pulse sequence and plotting the PL inten¬ 
sity as a function of trf- Typical nuclear spin Rabi oscil¬ 
lations with B rf = 12.5 G are shown in Fig. [2jd) . The 
nuclear Rabi frequency fljv = 48 kHz we extract from 
the data exceeds the expected value = gN^NBpF = 
3.85 kHz by more than a factor of 10, indicating an en¬ 
hancement of the effective nuclear gyromagnetic ratio. 

We further examine the enhancement by plotting the 
nuclear Rabi frequency Ojv as a function of B rf for a 
series of dc magnetic fields B z (see Fig. 3). For each 
value of B z , Qjy scales linearly with Hr p, with the ob¬ 
served slope implying a constant effective gyromagnetic 
ratio 7jv, e ff whose value is greater than the bare nuclear 
gyromagnetic ratio yjv = gNUN = 0.308 kHz/G. The 
data also confirms the dc magnetic field is well aligned 
with the NV axis, as there is no visible offset at Rrf = 0 
that would result from an off-axis magnetic field B x y . 
The enhancement factor ")N,es/lN = 18 at B z = 514 G 
and increases to 57 at B z = 873 G. 

The enhancement of the effective nuclear gyromagnetic 
ratio can be understood by examining the secular and 
non-secular terms in the Hamiltonian, which are given 


E^Jh = - 


g e H B Aj 


D - g e \i B B 
The nuclear Rabi frequency is then given by: 


( B X I X + Byly). (6) 


Hat = 


(gN^N - 


g e y B A ± \ 

D - g e ^ B B z ) 


-Brf, 


( 7 ) 


where Hrf is perpendicular to the NV axis. 

It is now apparent that second order perturbation the¬ 
ory accounts for the enhanced effective nuclear gyromag¬ 
netic ratio observed in the experiment: 


TAT.eff = gNUN — 


g e H B Aj_ 

D - g e ^ B B z ' 


( 8 ) 


The enhancement of the nuclear gyromagnetic ra¬ 
tio is analogous to two-qubit interactions in the cir¬ 
cuit quantum electrodynamics (cQED) architecture 
[Fig. [4l(a) ] 27 i 28 . In cQED, two qubits can be coupled via 
the exchange of virtual photons in a microwave cavity. 
The effective coupling rate J ~ c/i^/AcQed, is a func¬ 
tion of the individual coupling strength of each qubit to 
the cavity (described by g\ and < 72 ) and the qubit-cavity 
detuning A c qed^. Here our NV electronic spin acts 
in place of the cavity, mediating the coupling between 
the nuclear spin (qubit 1, coupling strength g\ ~ A±) 
and the RF magnetic field (“qubit” 2, coupling strength 
52 ~ gefi B )■ The detuning is set by the energy difference 
A between the ms = 0 and ms = — 1 states. 

In Fig. 4(b) we plot 7iv,eff as a function of 1/A. The 
data are extracted from Fig. 3 and the enhancement of 
the effective nuclear gyromagnetic ratio is in good agree¬ 
ment with perturbation theory (solid line). We observe 
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FIG. 4: (a) Comparison of the mechanisms that lead to the 
cavity-mediated couplings of two qubits in the cQED archi¬ 
tecture and the nuclear spin gyromagnetic ratio enhancement 
in NV centers (see text for details), (b) Effective nuclear spin 
gyromagnetic ratio, 7;v,eff, plotted as a function of 1/A, and 
extracted from the nuclear spin Rabi data shown in Fig. [3] 
Enhancements 7jv,eff/7JV > 50 are obtained. The data are 
in good agreement with the theoretical value given by Eq. [8] 
(solid line). 


enhancements 7jv,eff/7iv exceeding 50, making it possi¬ 
ble to perform a full nuclear spin Rabi oscillation on a 
few microsecond timescale with only ~ 10 mW of power 
delivered to the sample. This enhancement is approxi¬ 
mately 5 times greater than earlier results obtained on 
NVs at lower magnetic fields^. We expect the observed 
trends to hold for smaller energy splittings than are used 
in this experiment, with the trade-off being the reduc¬ 
tion in the nuclear spin lifetime T) jv near the ground 
state level anti-crossing^, where the energy splitting ap¬ 
proaches zero. These smaller energy splittings are not 
accessible with our current experimental setup due to 
the difficulty of polarizing the nuclear spin far from the 
ESLAC. 

In conclusion, we have shown that the effective nuclear 
spin gyromagnetic ratio can be greatly enhanced due to 
hyperfme coupling of the nuclear spin to the NV elec¬ 
tronic spin. The enhancement is well described by sec¬ 
ond order perturbation theory and is analogous to cavity- 
mediated qubit couplings in cQED. Our approach is also 
applicable to other coupled electron-nuclear spin systems, 
such as phosphorous donors in silicon or rare-earth metal 
dopants in yttrium orthosilicat o 10 ’ 12 , where it would al¬ 
low for rapid quantum control of nuclear spins without 
requiring high RF power. 
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